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THE  VIRTUAL  MASS  AND  LIFT  FORCE  ON  AN  OBLATE 
ELLIPSOID  IN  ROTATING  AND  STRAINING  INVISCID  FLOW 
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Abstract:  The  introduction  of  a  rotating  frame  allows  the  use  of  a  potential  fiinction  to 
solve  the  problem  of  a  ri^d  oblate  ellipsoid  undergoing  rotation  and  translation  in  a 
rotating  and  straining  inviscid  flow.  Virtual  mass  and  lift  force  coefficients  are  calculated, 
including  a  coefficient  associated  with  the  ellipsoid  rotation  contribution  to  the  lift  force. 
Significantly,  the  computed  coefficients  can  depart  from  sphere  values. 
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1.  Introduction. 

The  virtual  mass  and  the  lift  force  on  a  sphere  has  been  derived  by  several  authors 
using  inviscid  flow  theory.  The  original  calculation  of  the  so-called  virtual  mass  force  for  a 
sphere  is  attributed  to  Lord  Kelvin  [Lamb,  1932].  This  force  arises  when  a  rigid  body 
accelerates  through  a  quiescent  fluid. 
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The  lift  force  i.e.,  the  force  on  a  body  moving  relative  to  a  straining  fluid,  has  been 
calculated  by  Auton  (1987)  and  Drew  and  Lahq^  (1987)  for  the  particular  case  of  a  rigid 
sphere.  In  order  to  make  the  problem  linear,  Auton  assumed  that  the  velocity  change  on 
the  sphere’s  surface  due  to  the  rotation  of  the  fluid  was  much  smaller  than  the  relative 
translation  velocity  of  the  sphere.  Drew  and  Lahey  (1987)  introduced  a  rotating  firame  of 
reference  to  eliminate  the  vorticity  and  to  be  able  to  use  a  scalar  potential  function  to 
calculate  the  velocity  field.  Since  the  vorticity  is  a  time-dependent  magnitude,  their 
solution  is  valid  for  the  interval  for  which  the  vorticity  in  the  rotating  fi'ame  is  small 
compared  with  the  fluid  vorticity  at  infinity  in  the  inertial  non-rotating  system.  Drew  and 
Lahey  (1990)  obtained  the  same  result  than  Auton’ s  and  demonstrated  that  their  approach 
is  more  general. 

When  a  gas  bubble  moves  in  a  straining  liquid  it  may  undergo  deformatioa  Thus, 
even  if  the  bubble  was  initially  a  sphere,  it  will  change  its  form  to  one  that  makes  all  the 
calculations  mentioned  above  non-applicable,  since  the  bubble  is  no  longer  a  sphere  nor  a 
rigid  body.  Moreover,  there  is  evidence  that  the  failure  to  model  bubbly  flow  after  a 
sudden  expansion  is  due  to  the  fact  that  lift  force  modeling  considers  only  rigid  spherical 
bubbles  [Bel  Fdhila,  1991].  Spite  its  practical  importance,  very  little  attention  has  been 
given  to  the  force  over  non-spherical  bodies.  The  authors  are  aware  of  only  the  work  of 
Lamb  (1932),  who  calculated  the  virtual  mass  coefBcient  for  an  ellipsoid  moving  in  a 
quiescent  fluid. 

It  is  the  purpose  of  this  report  to  calculate  the  virtual  mass  and  lift  force  on  a  rigid  - 
oblate  ellipsoid  (i.e.,  a  ellipsoid  obtained  by  the  rotation  of  an  ellipse  around  its  minor  'f 
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semiaxis).  Drew  and  Lahey’s  approach  is  used  in  section  2.  In  section  3  the  velocity 
potential  is  calculated.  Finally,  in  section  4  the  force  on  an  ellipsoid  is  obtained  and  from 
it,  the  appropriate  coefScients  are  derived. 

2.  Statement  of  the  problem  and  theoretical  approach. 

The  problem  to  be  solved  is  shown  schematically  in  Figure- 1.  An  oblate  ellipsoid, 
whose  center  is,  using  indicial  notation,  at  position  x*j,  moves  through  an  inviscid  fluid 
with  an  unsteady  translation  velocity  of  v*ji  and  rotates  at  a  constant  angular  velocity  of 
a)*di  (the  superscript  *  indicates  the  inertial  system,  the  subscript  d  stands  for  the  discrete 
phase  and  the  dummy  index  i  indicates  the  Cartesian  components).  The  fluid  far  from  the 
ellipsoid  is  undergoing  a  motion  which  consists  of  unsteady  translation,  a  rotation  at 
constant  angular  velocity  and  a  constant  str^.  The  velocity  of  the  continuous  phase  far 
from  the  ellipsoid  is  thus 

V  *ci  =  V  *0i  %•  JC  *i  +e,^0  *o;  X  **  (1) 


where  x*i  is  the  spatial  coordinate,  v*o,  is  the  time-dependent  undisturbed  fluid  velocity 
(the  subscript  0  stands  for  undisturbed),  e%=  (v'''e<^+  v*g;,)/2  is  the  constant  symmetric 
tensor  representing  the  straining  motion  and  (O*oi  is  the  constant  angular  rotation  vector. 
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In  order  to  get  the  velocity  field  fi'om  just  a  scalar  velocity  potential  it  is  necessary 
to  eliminate  the  vorticity.  For  that  purpose  a  noninertial  (unstarred)  coordinate  system 
which  rotates  with  the  fluid  far  fi’om  the  ellipsoid  is  introduced.  The  systems  are  related  by 


where  Qij  is  an  orthonormal  rotation  tensor  which  satisfies 


(2) 


QijQig  -Qji  -  ^ik 


(3) 


and 


(4) 


It  can  be  concluded  fi’om  Eq.  (2)  that 


Vi  =  QijV*j+Qij  x*j  (5) 

For  an  incompressible  inviscid  fluid  in  a  gravitational  field  a  vorticity  balance  gives 
[Drew  and  Lahey,  1990] 


C>,+CmV*-(2a  +  C*K>=0 


(6) 
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where  Q  -Ejik  v*,  is  the  vorticity  in  the  rotating  unstarred  system.  From  Eq.  (6)  it  can  be 
seen  that  the  vorticity  is  not  a  constant  magnitude,  but  depends  on  time.  Thus  the  new 
coordinate  system  can  have  zero  vorticity  at  a  given  instant,  chosen  arbitrarily  to  be  ^=0, 
but  later  it  will  grow.  As  a  consequence  the  present  analysis  is  valid  only  in  the  time  scale 
that  maintmns  the  vorticity  at  a  small  enough  value  for  the  potential  theory  approximation 
to  be  satisfactory  [Drew  and  Lahey,  1990]. 

In  order  to  simplify  the  force  calculation  it  is  convenient  to  derive  the  potential  in  a 
non-rotating  system.  Thus,  the  following  change  of  variables  is  introduced; 

Xi  =  Qa  Xj  (7) 

=  Qji  V,  (8) 

Notice  that  in  the  new  (hat)  system,  the  velocities  are  not  obtained  as  the 
derivatives  of  the  position.  The  hat  and  star  system  are  related  by 

x,=x*,  (9) 

Vi  =v*i  +QjiQj^x*^ 


(10) 


Replacing  Eq.  (1)  in  Eq.  (10)  and  using  Eq.  (4),  we  may  conclude  that: 


=  v*„;  +e*>.  X*; 


(11) 


Thus  the  motion  of  the  fluid  far  from  the  ellipsoid  consists  of  an  unsteady  translation,  at 
velocity  Voi  =  and  a  constant  strain  -  The  velocity  on  the  ellipsoid’s  surface 

can  be  obt^ed  in  a  oimilar  manner; 


Vv;  =  V  ’ 


+e.a 


[(a%+GJ%l 


(12) 


where,  in  the  last  equation,  belongs  to  the  ellipsoid’s  surface  and  (0*ef  is  the  ellipsoid 
rotation  in  the  inertial  (starred)  system.  From  the  previous  equations  it  is  clear  that  an 
observer  attached  to  the  hat  system  see  the  ellipsoid  undergoing  an  unsteady 
translation  at  a  velocity  v*^,.  and  a  rotation  with  an  angular  velocity  This 

rotation  means  that  the  ellipsoid’s  surface  is  moving.  This  movement  will  make  the 
boundary  conditions  on  the  ellipsoid’s  surface  dependent  on  the  surfrice  position,  except 
for  the  particular  case  in  which  0}%+(o%  is  parallel  to  the  direction  of  the  minor  semiaxis 
of  the  ellipsoid.  It  will  be  seen  later  that  there  is  an  analytical  solution  for  the  velocity 
potential  that  satisfies  the  rotation  boundary  condition. 


3.  The  velocity  potential  in  the  hat  system. 


The  velocity  potential,  (f>,  is  defined  by, 

V,  =  0,  (13) 

and  consequently,  for  an  irrotational  flow  satisfies  the  Laplace  equation, 

0^=0  (14) 

Oblate  spheroidal  coordinates  are  defined  by  [Morse  and  Feshbach,  19S3]; 

Jci  -  icji  =  a  COS0  (15) 

Jcj  - jCj2  =  a  +1  sin© 

where  a>0  is  a  constant,  §  goes  firom  0  to  infinity,  q  goes  fi’om  -1  to  1  and  fl  goes  firom  -n 
to  7t.  All  the  magnitudes  and  operators  required  to  calculate  the  force  on  the  ellipsoid  in 
the  oblate  system  of  coordinates  can  be  derived  firom  these  equations.  The  more  important 
magnitudes  are  presented  in  the  Appencfix.  The  surfaces  |  =  >  0  are  oblate  ellipsoids 

with  major  and  minor  semiaxes  given  by. 
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b  =  a  ^l  +  lo^  and  c  =  a  |o, 


(16) 

respectively.  From  the  previous  equations  it  can  be  seen  that  a  sphere  of  radius  R  is  the 
limiting  case  obtained  when  goes  to  infinity  while  the  product  R  =a^o  is  kept  constant. 
The  surfaces  t)  =  Ho  are  hyperboloids  of  one  sheet,  asymptotic  to  the  cone  of  angle 
arcos(T])  with  respect  to  the  Jcj  axis,  which  is  the  axis  of  the  cone. 

The  general  solution  for  the  Laplace  eaquation  in  this  coordinates  system  is  [Morse 
andFeshbach,  1953]: 

♦(? .n .«)  =  co^me )  +  sm(m9 )] P,” (17) [C„  C (H )*D„a: (i| )| 

IVsO  JTPsO 

(17) 

where  P”  and  2"  are  the  Legendre  fimctions  of  first  and  second  kind  respectively, 
A^,  are  coeflBcients  that  depend  only  on  indices  m  and  n  and,  i  = 

It  must  be  noted  that  the  coefiBcients  Cmn  and  may  be  complex  numbers  as  necessary 
to  make  the  velocity  potential  a  real  fimction.  To  avoid  this  unnecessary  complication,  the 
Legendre  fimctions  used  in  this  work  were  redefined  in  such  a  way  that  only  real  constant 
coefficients  were  needed.  These  Legendre  fimctions  are  listed  in  the  Appendix . 

In  order  to  calculate  the  potential,  four  cases  are  considered:  In  the  first,  the 
ellipsoid  4  =  >  0  is  moving  with  velocity  v*ji  and  the  fluid  at  infinity  is  at  rest;  in  the 
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second,  the  ellipsoid  is  at  rest  and  the  fluid  far  away  fi-om  the  ellipsoid  is  moving  with  a 
velocity  v*oi',  in  the  third  case,  the  velocity  at  infinity  consists  of  the  strain  tensor  and  the 
ellipsoid  is  at  rest;  and  in  the  last  case,  the  ellipsoid  is  undergoing  a  rotation  and  the  fluid 
at  infinity  is  at  rest.  The  general  solution  is  formed  by  superposition  of  these  four 
particular  ones. 

The  boundary  conditions  for  the  first  case  are,  on  the  ellipsoid’s  surface; 

08) 

where  rij  is  the  normal  to  the  ellipsoid’s  surface,  and  at  infinity: 

<t>  -j^O  (19) 

In  order  to  find  the  velocity  potential,  the  velocity  v*ji  is  separated  into  its  Cartesian 
components  and  the  gradient  is  expressed  in  the  oblate  coordinate  systems  to  obtain 
[Lamb,  1932]; 


0  =  a[v*^i  COS0 


+av*,j 


(20) 
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The  boundary  condition  for  the  second  case,  at  injSnity  is; 


and  at  the  ellipsoid’s  sur&ce; 


(21) 


<t>.jnj=0  (22) 

The  solution  can  be  found  by  replacing  Eqs.  (15)  in  Eq.  (21)  and  comparing  the 
result  withEq.  (17).  Since  the  functions  Q"  (i§)  tend  asymptotically  to  zero  when  ^  tends 
to  inifinity,  the  functions  P”  (i^ )  are  the  ones  that  satisfy  the  boundary  condition  at  infinity 
and  the  functions  Q"  (i§)  are  useful  for  the  ellipsoid’s  surface  boundary  condition.  Thus, 
the  solution  is: 
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» 


^  =  a[v*oi  COS0  +v*o,  sine] 


?•+ 


+aVo3Pi°(n) 


(23) 


The  boundaiy  conditions  for  the  third  case  are,  at  infinity; 


(24) 


and  at  the  ellipsoid’s  surface  the  boundaiy  condition  is  the  same  as  that  for  the  second 
case,  Eq.  (22).  The  same  procedure  that  was  already  used  for  the  second  case  can  be 
applied  here  to  demonstrate  that  a  solution  that  involves  only  a  finite  number  of  terms  in 
Eq.  (17)  is  possible  only  for  the  particular  str^  tensor  given  by. 
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^  *23  » 

0 


(25) 


and  the  solution  is; 


The  boundary  condition  foi  the  fourth  case  at  the  ellipsoid’s  surface  comes  from 
the  second  term  of  Eq.  (12)  and  describes  a  rotation: 


<pjnj=e^{(u%HO%}x,n, 


(27) 
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The  condition  at  infinity  is  the  same  as  for  the  first  case  (see  Eq.  (19)).  Assuming  that  the 
ellipsoid’s  minor  semiaxis  is  parallel  to  the  Jc,  axis,  (0*(^(0*e^  can  be  expressed  in  terms  of 

the  oblate  unit  vectors  and  the  right  hand  »de  of  (27)  expressed  in  terms  of  the  oblate 
coordinate  system.  It  would  seem  that  the  expression  obtained  in  this  way  is  valid  only  at 
the  instant  for  which  the  ellipsoid  minor  semiaxis  is  parallel  to  the  axis,  because  later 
the  ellipsoid’s  rotation  will  change  the  expression  of  Q)*(j+Q)%  in  terms  of  the  oblate  unit 
vectors.  However,  since  0)*of+(0*^  is  parallel  to  the  axis  of  rotation  of  the  ellipsoid,  its 
projection  on  the  oblate  unit  vectors  at  a  given  point  over  the  surface  does  not  change 
after  a  rotation.  Thus,  the  expression  obtained  assuming  a  particular  orientation  of  the 
ellipsoid  with  respect  to  co*(^+ct}*e^  is  completely  general.  The  result  is: 


-H 


CO  *01  +0)  * 


\n)sme  +(i 


)  COS  8]  77^1^ 


^|^o+r 


(28) 


As  expected,  the  co*oi+co*di  component  does  not  appear  in  equation  (28).  Since  the  Jcj 
axis  is  a  rotational  symmetry  axis,  when  the  ellipsoid  rotates  around  the  x,  axis  it  does  not 
displace  liquid.  Thus,  for  this  particular  rotation  the  boundary  condition  is  similar  to  that 
of  the  sphere  which  is  independent  of  its  rotation.  Combining  Eqs.  (28)  and  (27),  and 
expressing  the  left  hand  side  of  this  last  equation  in  the  oblate  coordinate  system,  the 
solution  for  this  case  is  obtained  as. 
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(29) 


Interestingly,  Lamb  first  derived  this  solution  in  1932. 

The  final  solution  for  the  problem  is  obtained  by  adding  the  four  particular  ones. 
However,  since  the  solutions  for  the  first  three  cases  assume  an  orientation  of  the  ellipsoid 
with  respect  to  the  inertial  system,  in  order  to  apply  the  superposition  principle  to  build  a 
general  solution,  v*oi,  v**  and  e*^  have  to  be  projected  to  a  Same  of  reference  attached  to 
the  ellipsoid.  This  is  accomplished  by  replacing  v*oi^  v*di  and  e*g  in  the  already  derived 
solutions  for  the  four  particular  cases  (i.e.,  Eqs.  (20),(23),(26)  and  (29))  by. 


(30) 


(31) 


and. 


(32) 


respectively,  where. 


with. 


4('=»)=«s 


A 

In  order  to  simplify  the  notation,  from  now  on  Q^-  will  be  dropped  from  the 

equations.  Thus,  the  forces  to  be  calculated  in  the  next  section  will  be  valid  at  the  initial 
instant  /=0.  To  obtain  the  forces  at  any  moment  later,  v*o,,  v*^,  and  e%will  have  to  be 
replaced  by  their  values  given  in  Eqs.  (30)  to  (32).  With  all  these  qualifications,  the  final 
velocity  potential  can  be  written  as. 


^  =  a[v*j,cos9  +  v%j  sine]  ;;'(n) 


«o  e;(i«) 


■+av*j3Pi“(7]) 


Gf(«'§) 


■‘u 


Ld| 


+a[(v*oi  +JC*«3  «*i3)cos0  +(v*02  +j:*^3  *23)0110]  P/(t]) 


L|e;0«) 


e!(ii) 


+a(v*03  e*i3  +.«*i2 «  *23)^1°  (n)i 


* 

A 

a“(ii)  • 

+Tt®*i3  +e*23Sin0]P2(Tj) 


P2W- 


Glte) 


-•{o 


I  + 


+y[-(QJ*oi  ■«w*Jsin0  +(g)*02  +®  %2)cos0|P2‘(t})-j|^ 


Mii) 


Id^ 


(33) 


4.  Derivation  of  the  force  on  an  ellipsoid. 


The  momentum  balance  in  the  inertial  system  is  [Morse  and  Feshbach,  1953]: 
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where  p  is  the  continuous  phase  denaty. 


Integration  of  the  previous  equation  along  a  current  line  from  infinity  to  the 
ellipsoid’s  surface  ^  yields: 


_(p  (”)  P  +  |[(v*,  v*.)(oo)-(v  V  (35) 


where  dl*j  is  the  arc  length  differentials  along  the  current  line.  It  can  also  be  seen  that  the 
integral  of  the  last  term  in  Eq.  (34)  drops  out  because  £^v  *jdl\  =  0. 

In  order  to  evaluate  the  right  hand  side  of  Eq.  (35),  it  is  convenient  to  express  the 
velocity  of  the  inertial  system  as  a  function  of  that  in  the  hat  systeuL  In  particular,  from 
Eq.  (10)  and  (4),  and  the  fricts  that  Q}*j  is  constant  and  the  time  derivative  is  partial,  it  is 
concluded  that, 


The  velocity  potential  was  derived  with  the  ellipsoid  at  the  origin.  Thus,  the  velocity 
potential  is  a  function,/,  of  Jcj  -  and  t : 


(37) 


v^ere  the  potential  depends  explicitly  on  time  through  and  v^.  ft  is  concluded  from 
Eqs.  (13),  (36).  (37)  and  (9)  that. 


(38) 


From  this  last  relationship  it  is  readily  seen  that. 


Iv*„  rfi*  =(#,  -f;V  %)(«)-(«,  (39) 

it 


After  replacement  of  Eq.  (39)  in  Eq.  (35)  and  rearrangement,  the  following  relationship  is 
obtained: 


where  the  pressure  at  infinity  was  chosen  in  such  a  way  that  Eq.  (40)  holds. 

Before  integration  of  the  last  equation  over  the  ellipsoid’s  surface,  h  is  convenient 
to  use  Eqs.  (10)  and  (4)  to  obtain: 


v*jv*j  =  VjVj+2eju<o *0*  X*,  v^.  \  x*,  (0*^x\ 


(41) 


IS 


The  net  force  on  the  ellipsoid  is  then  obtained  as  the  integral  of  Eq.  (40)  over  the 
ellipsoid’s  sur&ce; 


-iip*nidS=  pii<p/iidS- pj\<pjV\jnidS+-jilv*jV*jnidS  (42) 

It  should  be  noted  that  rii  is  the  unit  normal  at  the  ellipsoid’s  surface.  The  area  differential 
dS  can  be  found  in  the  Appendix.  The  last  integral  in  Eq.  (42)  can  be  expanded  using  Eq. 
(41): 

I  JJv  * .  V  =  I  If  VjVj  n^dS+\\e^G)  (x  *,  -x \)vj  dS  (43) 


The  integral  over  the  ellipsoid’s  surface  of  the  last  term  in  the  right  hand  side  of  Eq.  (41) 
is  zero  because 

JJ(x\-xV)(jf*.-Jc*JnidS  =  0,  (44) 

as  can  be  demonstrated  multipl^g  the  integral  by  a  constant  vector  and  iq)piying  the 
divergence  theorem  to  the  result.  This  result  is  an  extension  of  a  similar  one  used  for  a 
sphere  [Voinov,  1973]. 
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The  following  result  for  the  net  force  on  the  ellipsoid  is  obtained  by  combining 


Eqs.  (43)  and  (42): 


-If P  *  ^  =  pJJ ^  -  pif ^ j V  %  Hi  iiS  +  y If VjVj  Tii  dS 

+  pUbjuO) *ot  (jr*,  -x*a)vj n,  dS 


(45) 


Using  the  velocity  potential  function,  0,  calculated  in  the  previous  section,  all  the 
integrals  in  the  last  equation  can  be  evaluated.  To  make  the  exposition  of  results  simpler, 
from  now  on  it  will  be  assumed  that  xjt  =  0.  To  obtain  the  force  for  the  case  xjt  ^  0,  the 
term  jc  has  to  be  added  to  v  V  The  results  are: 


V 


* 


(47) 
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(50) 


.  /3  =  §0  (lo"  + 1)  J-  ^0* 

It  is  worth  noting  that  the  integral  in  Eq.  (48)  is  the  only  one  where  the  ellipsoid 
rotation,  (o*dt,  provides  a  non-zero  contribution  to  the  force.  From  Eqs.  (45)  to  (49),  the 
net  force  on  the  ellipsoid  can  be  written  as, 

-iip*rt,clS=  p  +l)  + 

~  2  -  *01  *<a)  ^  *0k  *03  *d%) 

®  %3  -  V  *03)_  + 

«•  -1 

+  P  +l)j[v*a.r  -«**  V**  +e\v\,  -ea,©*0*  V%,]  (/  =  1,2) 

=  P  +l)]  ^(v*i«.,-vVr)+^^«*«(  V  V  -  V**) 

~  2-  {  ''*ot  2-  ®%(  ''*<*  ~''*ok) 

+ P  (§0  + 1)]  [v  *oi.t  -e  *»  V  •*  +e  V  *0*  -e^co  v  *„]  (/  =  3) 

(51) 

It  is  convenient  to  note  that. 
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where />*o  is  the  pressure  in  the  undisturbed  fluid.  The  first  two  terms  in  Eq.  (51)  are  the 
virtual  mass  force, 


(/  =  3) 
(53) 


Lamb  (1932)  obtained  this  result  for  the  particular  case  Since  P  tends  to 

2/3  when  goes  to  infinity,  the  well  known  coefficient  value  of  1/2  [Auton,  1987]  [Drew 
and  Lahey,  1987]  is  recovered  for  the  sphere  limit.  The  coefficients  P  f{2~  fi) 
and(l  -  P)/P  can  be  seen  in  Figure-  2.  To  interpret  this  figure  it  is  convenient  to  consider 
the  following  relationship,  which  can  be  derived  firom  Eqs.  (16), 


(54) 


As  expected,  for  the  special  case  of  e**  sO,  the  virtual  mass  coefiBcient  in  the 
direction  is  always  greater  than  for  the  directions  /=1,2.  This  results  because  the  area 
projected  by  the  ellipsoid  into  a  plane  perpendicular  to  the  direction  of  movement  is 
proportional  to  for  i=2  and  to  bc<  b^  for  /=1,2.  Thus,  the  amoimt  of  fluid  displaced  in 
the  former  case  has  to  be  greater  than  in  the  latter.  Since  for  both  directions,  /=1,2,  the 
virtual  mass  co^cient  is  the  same  (when  ^**=0),  the  problem  of  obtaining  the  virtual 
mass  coefiBcient  for  an  arbitrary  direction  is  bidimensional.  Thus,  the  coefiBcient  of  virtual 

volume,  Cmi,  when  the  relative  acceleration  of  the  ellipsoid,  forms  an 

angle  9)  with  the  axis  is. 


A  similar  definition  is  obtained  for  the  case  (v  ,  -v  ♦  s  0  and  e*^^. 

Even  though  the  virtual  mass  coefiBcient  (l  -  )/j9  tends  to  infinity  when  ^0  tends 

to  zero,  the  total  virtual  mass  force  is  not  infinite  because,  if  a  is  kept  constant,  the 
ellipsoid  volume  tends  to  zero.  The  total  force,  which  corresponds  to  a  fiat  disk  of  radius 
a,  can  be  obtained  firom  the  following  limit. 
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When  the  ellipsoid  volume  is  kept  constant,  the  total  force  and  a  go  to  infinity 
when  |o  tends  to  zero.  Since  a  is  the  disk  radius,  the  total  force  needed  to  accelerate  the 
disk  has  to  diverge. 

The  last  two  terms  in  Eq.  (S 1)  are  the  lift  force, 

(/  =  1,2) 

(/  =  3) 
(57) 

Taking  the  limit  when  tends  to  2/3,  it  can  readily  be  shown  that  the  sphere’s  lift 
force  is  recovered  from  the  last  equation.  There  are  two  different  contributions  in  the  last 
equation,  one  coming  from  the  vorticity  &r  away  from  the  ellipsoid,  (O*oh  and  the  other 
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coining  from  the  ellipsoid’s  rotation,  ai*di-  Let  us  consider  the  former  contribution  first; 
When  the  relative  velocity  -{v  *o,.  is  in  the  direction  /=3  (or  equivalently  (0*q/is  in 
the  direction y=l,2),  the  lift  force  coefficient  is  (l-  ^)//3  .  Otherwise  the  coefficient  is 
/(2  -  ^ .  It  must  be  noted  that,  due  to  the  same  reason  as  in  the  virtual  mass  case,  the 
bigger  lift  force  coefficient,  (l-  ^)/j3  ,  occurs  when  the  projected  area  of  the  ellipsoid 
into  a  plane  perpendicular  to  the  relative  velocity  is  the  greater  and  the  smaller  lift  force 
coefficient,  ^ /(2 -  ^),  occurs  in  the  opposite  situation.  The  problem  of  determining  the 
lift  force  coefficient  for  an  ellipsoid  moving  with  a  relative  velocity,  -{v  *0. in  an 
arbitrary  direction,  is  also  bidimensional  and  the  coefficient  is  the  same  as  in  Eq.  (55), 
where  now  <p  is  the  angle  between  the  relative  velocity  and  the  axis.  It  is  also  clear 
that,  in  comparison  with  the  sphere’s  lift  force,  neither  forces  in  other  directions  nor 
changes  in  sign  arise.  Since  the  1/2  value  for  the  sphere’s  case  is  reached  asymptotically  for 
both  coefficients  (see  Figure-  2),  the  use  of  the  sphere  value  is  acceptable  in  any  situation 
where  bubbles  are  not  subject  to  strong  deformation  forces.  However,  if  departure  from 
the  spherical  shape  is  important,  for  example  when  bubbles  are  subject  to  a  strong  shear, 
or  acceleration  is  not  negligible,  the  shape  effect  becomes  important.  For  instance,  a 
bubble  with  b=2c  has  ^=1/3^,  which  implies  ^  /(2  -  ^ )  w  0.31  and  (l -  ^ )/j3  »  U5, 
and  the  predictions  are  off  apprommately  by  a  frctor  2  ,  from  that  of  a  sphere. 

Figure-  2.  As  for  the  (O*o,  terms,  the  coefficient  corresponding  to  the  relative  velocity  in 


can  also  be  seen  in 


The  coefficients 


f'A-yi 


and  : 
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the  direction  i=3  is  the  bigger.  As  expected,  the  component  (0*j3  does  not  contribute  to 
the  lift  force.  It  must  be  noticed  that  depending  upon  the  relative  values  of  co**,  and 
/3  the  dominating  terms  might  be  the  ones  originated  by  the  ellipsoid  rotation,  (o*dt.  These 
terms  could  be  important  in  geometries  of  practical  interest,  such  as  a  sudden  expansion 
where  viscous  efiects  will  make  the  deformed  bubble  rotate. 

5.  Conclusions. 

The  virtual  mass  and  the  lift  force  for  a  rotating  oblate  ellipsoid  moving  in  a 
in^dscid  fluid  undergoing  an  unsteady  translation,  constant  strain  and  constant  rotation, 
were  calculated.  The  lift  force  includes  the  contribution  of  the  ellipsoid’s  rotation.  The 
coeflBcients  of  virtual  mass  and  lift  force  differ  appreciably  fi-om  those  of  the  sphere  only 
for  ellipsoids  of  considerable  eccentricity.  Thus,  it  is  acceptable  to  use  the  sphere 
coefiicients  when  modeling  two-phase  flow  regimes  with  no  important  deformation 
effects.  In  contrast,  whenever  there  are  significant  deformation  effects  in  bubbles  shape,  it 
is  more  appropriate  to  use  the  coefidcients  obtained  in  this  work. 

One  of  the  more  striking  results  is  that  the  virtual  mass  and  lift  force  coeflBcients 
for  a  given  relative  velocity  direction  are  the  same,  although  there  is  no  evident  reason  for 
this  to  occur.  It  seems  worth  it  to  investigate  this  issue  further. 
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6.  Appendix:  Properties  of  oblate  coordinates  and  Legendre 
functions  definitions. 


Unless  explicitly  stated,  all  definitions  were  taken  fi-om  the  woric  of  Morse  and 
Feshbach  (1953). 


Scale  factors: 


Gradient: 


where  I",  if  and  6  are  the  unit  vectors  normal  to  the  surfeces§  =  |o*  H  =^0  6  =  6^, 

respectively,  where  |o,  rjo  and  do  are  constants. 

Area  differential  over  the  ellipsoid’s  sur&ce: 


dS  =  h^\de  dr\ 


Legendre  Fimctions: 
Functions  of  a  real  argument; 


PiM  =  3ti 

Functions  of  an  imaginary  argument,  as  redefined  in  this  work  to  avoid  imaginary  numbers; 

p, "(<«)=« 

?,'(<«  )=Vi+?' 

?,'(/«)= 35 


a"(i|)  =  5  aretan(>^)-l 

a'('0=  -Vl+S’ 


7.  References. 


Auton,  T.  R..  “The  lift  force  on  a  spherical  body  in  a  rotational  flow”.  J.  Fluid  Mech.  183, 
pp.  199-218, 1987. 

Bel  Fdhila,  “Analyse  Experimentale  et  Modelisation  d’un  6coulement  vertical  a  buUes  dans 
un  ilar^ssement  brusque”.  PhD  Thesis,  Universite  de  Toulouse,  Toulouse,  France,  1991. 

Drew,  D.  A.  and  Lahey,  R.  T.  Jr.  “The  virtual  mass  and  lift  force  on  a  sphere  in  rotating 
and  straining  inviscid  flow”.  Int.  J.  Multiphase  Flow.  13,  No.  1,  pp.  113-121,  1987. 

Drew,  D.  A.  and  Lahey,  R.  T.  Jr.  “Some  supplemental  analysis  concerning  the  virtual  mass 
and  lift  force  on  a  sphere  in  a  rotating  and  straining  flow”.  Int,  J.  Multiphase  Flow.  16, 
No.  6,  pp.  1127-1130,  1990. 

Lamb,  H.  “Hydrodynaimcs”  Dover  publications,  New  Yoric,  1932. 


Morse,  P.  M.  and  Feshbach,  R  “Methods  of  Theoretical  Physics”  McGraw-Hill,  1953. 


Voinov,  0.  V.  “Force  actmg  on  a  sphere  in  an  inhomogeneous  flow  of  an  ideal 
incompressible  fluid”.  Zh.  prikl,  Melch.  tekhJSz.  14(4),  592-594. 


31 


Figare-l  :Description  of  the  problem  to  be  Mlved.  Vectors  or  tensors  are  bold. 

Kgure-  2:  The  coefficients  of  virtual  mass  and  lift  force  (see  Eqs.  (S3)  and  (57))  as  a  function  of 

where  b  and  c  are  the  major  and  minor  semiaxes  respectively  and 

(Eq.  (50)).  Notice  that  the  coefficients  that  tend  to  zero 
correspond  to  the  contribution  of  the  ellipsoid  rotation,  ,  to  the  lift  force. 
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figure-2 


